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I. INTRODUCTION 

The semiclassical theory of electronic dynamics has 
been successful in explaining a wide range of transport 
phenomena in solid-state physics. It is well established 
that the semiclassical single-particle equations of motion 
(sEOM) acquire anomalous corrections that stem from 
Berry phase^l accumulated in the adiabatic motion of 
a wave packet .I^^Si Such Berry-phase effects have been 
successful in explaining the anomalous Hall effect,!^ cor- 
rections to semiclassical quantization,'^ intrinsic m ag- 
netic moments of electronic wave p acket s,^ ^!^!^^'"'^^!^^ and 
anomalous thermoelectric transporlP^E^ various elec- 
tronic systems. The (noncanonical) Hamiltonian per- 
spective and issues related to the Liouville's theo- 
rem and modified electron density of states have also 
been elucidatedrecently,'SHIIl as well as Fermi-liquid 
generalizations.'^^ 

Collective dynamics in the semiclassical approach is 
usually described with the Boltzmann equation, where 
electrons drift in phase space according to their single- 
particle equations of motion projected on respective 
bands. While such an approach is intuitive and physically 
appealing, its validity is open to question. The problem 
is that spatiotemporal inhomogeneities in the Hamilto- 
nian such as electromagnetic potentials induce interband 
coherences that may not be completely captured by the 
sEOM. Furthermore, it was noted in Ref. TT] that the 
sEOM are noncanonical, appearing to violate Liouville's 
theorem. As a remedy, the authors introduced a rescal- 
ing of the phase-space density of states. This modifica- 
tion was shown to be consistent with a formal requan- 
tization procedure for band-projected Bloch electrons, 
which promotes the Poisson brackets of the sEOM to 
commutators, resulting in noncanonical commutator re- 
lations of position and momentum and the corresponding 
modified minimum quantum uncertainty in phase-space 
variables.!^ However, wc find no logical necessity to en- 
force Liouville's theorem for the band-projected distri- 
bution function, and the formal requantization argument 
does not explain why the density of states is modified by 
the band projection. 



These fundamental issues motivate us to consider the 
generalized Boltzmann equation for a multiband system 
from a systematic, ground-up approach. In this paper, 
for a clean system with nondegenerate bands, we derive 
a band-projected kinetic equation by performing a semi- 
classical expansion of the density-matrix equation of mo- 
tion. A similar Green's function approach was developed 
in Ref.ITSl However, there it was assumed from the out- 
set that the distribution function was nonvanishing only 
in a single band. In this paper, we make no a priori 
assumptions of decoupled bands and capture systemat- 
ically all corrections to the Boltzmann equation up to 
second order in space-time gradients. Despite our at- 
tempt at decoupling by a systematic gradient expansion, 
we find remaining interband terms that warrant further 
investigation. 



This paper is organized as follows. In Sec. |II A[ for 
a two-band Hamiltonian, we derive a covariant trans- 
port equation for the Wigner distribution function. In 
Sec. |II B{ wc decouple the transport equation to derive 
effective band-diagonal semiclassical kinetic equations 
(sKE). In our sKE, in addition to the known Berry curva- 
ture corrections to the sEOM, we find terms correspond- 
ing to the aforementioned modified density of states, 
as well as interband terms not seen before. We inter- 
pret these terms as representing quasiparticles motion 
in curved, nonorthogonal subspaces of the total Hilbert 
space. In Sec. |II C[ we extract the hydrodynamic cur- 
rent from the continuity equation, discuss issues associ- 
ated with momentum-space ultraviolet cutoffs and band- 
crossing points. In Sec. |IID] we consider minimal cou- 
pling to electromagnetic gauge fields and express our 
transport equations in terms of the gauge-invariant ki- 
netic momentum. We then apply these equations to 
the 2D Dirac and 3D Weyl Hamiltonians in Sec. |III A| 
and Sec. |III C) In the conclusion, Sec.|IV[ we qualita- 
tively compare our approach with other methods for de- 
riving band-projected effective Hamiltonians by canoni- 
cal transformations on the Hilbert space. Computational 
details are relegated to the appendices. 
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II. SEMICLASSICAL KINETIC EQUATIONS 
A. Covariant Transport Equation 

Consider the many-particle density matrix, 

Pap{ri,r2,t) = (V'|^(r2,0V'a(ri,i)) , (1) 

where (ri, t) are second-quantized field operators, and 
a is a spin or band index. We will consider only fermions 
in this paper, although our formalism is equally appli- 
cable to bosons. For a quadratic, mean-field (or Fermi- 
liquid) Hamiltonian, 

H{t)= J dridr2^piiri,t)nc.piri,r2,t)^,3ir2,t), (2) 

where H is the first-quantized, quasiparticle Hamiltonian 
expressed as a kernel (with the implied summation over 
repeated indices a,/?), the equation of motion for the 
density matrix closes: 

dtPal3 + ^{Ua-y P^I3 - Pa^ •S>'Hj0) = ■ (3) 

Here, ® denotes a convolution integral in real space: 
[A^ B](ri,r2) = J dr'A(ri, r')B(r', r2). We generally 
consider a smooth low-energy effective Hamiltonian 

■Ha/3(ri,r2,i) = ea/3(-i9ri , ri, i)(5(ri - . (4) 

If necessary, the function eQ^(— , ri, is properly 
symmetrized in its noncommuting arguments d^^ and ri . 
Eq. Q may represent the continuum Hamiltonian for the 
slowly-varying envelope fields in the k • p expansion, or 
the continuum limit of a tight-binding Hamiltonian. 

For slow and long-wavelength spatiotemporal inhomo- 
geneities, it is useful to define the distribution function 
np(r,t) by the Wigner transform (WT) of the density 
matrix: 

np(r,i)^ J dr'e'HP r'^(^r-f^,r-^,i) . (5) 

In the Wigner representation, the expectation value of 
an observable described by the kernel A(ri,r2) is given 
by 

(i)(t) - / dr^Tr[np(r,t)ip(r,t)] , (6) 
p 

where = / dp/ {2TThY, d being the number of spatial 
dimensions. Here and henceforth, we make a convention 
to denote the spin/band matrix structure by hats. When 
np(r, t) is a smooth function of r and t, one may construct 
a semiclassical kinetic equation for np{r,t) by taking the 
WT of Eq. ([3| and performing a gradient expansion. 

The kinetic equation for np(r,<) is governed by the 
quasiparticle energy matrix ep(r,t), which is the WT of 
the quasiparticle Hamiltonian kernel Q. Considering, 



for simplicity, a two-band system, the energy may be 
expanded as 

ep(r,<)- ^ep(r,t) + Ap(r,i).f . (7) 

Here, fa = (3'a/2 are the spin-1/2 matrices satisfying 
[Ta, h] = ie'^^'^Tc and {tq, ff,} = 5a6/2, (Ta being the Pauli 
matrices, e"*"^ is the antisymmetric Levi-Civita tensor, 
and we will use letters in the beginning of the alphabet for 
indices representing spin degrees of freedom. Below, we 
will be similarly decomposing the components of any 2x2 
matrix into scalar and vector pieces as M = M/2-|-M-f , 
where M = Tr[M] and M = Tr[Af(T]. We parametrize 
the gap vector as Ap(r,i) — Ap(r, t)mp(r, t), where 
Ap ~ I Ap I and nip = (sin 9^ cos ip-p , sin 6^ sin ip-p , cos 9^ ) 
is a unit-vector field represented by the spherical angles 
6'p(r,i) and i^p(r, i). 

The gap vector field Ap(r,t) appears formally as a 
magnetic field coupled to spin in phase space, whose di- 
rectional field mp(r,i) we will call the spin texture. Ex- 
amples of Hamiltonian ([?]) occur in ferromagnetic semi- 
conductors, where the gap vector represents the exchange 
and spin-orbit fields, and in nonmagnetic semiconductors 
with spatially varying spin-orbit coupling. The internal 
degrees of freedom need not be the actual electron spin. 
For example, our kinetic equation can be applied to the 
pseudospin dynamics near the K{K') points of graphene, 
which is described by the Dirac Hamiltonian. 

We will be interested in quasiparticle transport on the 
semiclassical spin-orbit bands defined by the eigenvalues 
of Eq.Q, 

ep.(r,t) = ^[ep(r,i) + sAp(r,0] • (8) 

where s = ±1. The local spin frame that diagonalizes 
Eq. I?! is defined by an SU(2) rotation C/p(r, t) such that 

f/t (Ap • T)C/p = Apf, . (9) 

The basis defined by this rotation is the local energy 
eigenstates in the sense that the average total energy {H) 
may be expressed as a spatial integral 

{H) = [ drY^^psnps- (10) 

ps 

Hps here are the diagonal elements of the distribution 
function in the local spin frame. 

To derive the semiclassical kinetic equation, we take 
the WT of Eq. ([3| up to second order in the gradient 
expansion. (See Appendix [A| for details.) The resulting 
kinetic equation in the local frame is 

i 1 
Dtfip + ^[ep,np] - -{Diip^D'hp} 

-'^[D,D,ep,D'D^np]^0, (11) 



3 



where {,} denotes anticommutator and [, ] commutator, 
i,j label coordinates of phase-space vector = (r,p), 
and summation over repeated indices is imphed. The in- 
dices are raised by 9* = J^-'dj, where J is the symplectic 
matrix acting on phase-space derivatives di = {dr, dp) as 



J 



J = 



1 

-1 



(12) 



the latter being written in terms of the d x d matrix 
blocks Pin Eq. (fTTl), we introduced the covariant deriva- 



tives of the distribution function and energy defined by 
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DiM = diM - i[Ai, M] = -diM + DiM 



(13) 



where ZJ/M 
tal letters /, 
phase-space 
gauge fields 
fined by Aj e 
where Af = 
the Euler-an; 



= djM. -f A/ X M. Hereafter, the capi- 
J, and K will be used to denote combined 
and time coordinates. The matrix-valued 
entering covariant derivatives 



(131 are de- 



-Ajr. 



Aj =F and f± = (f^ ± ify)/2. In 
gle parametrization of our local spin frame. 



e " e , the gauge fields are 



= cos edjip + dij , 



Af = -e^'''{smediip±idie), 



(14) 



where 7 is an arbitrary rotation angle about mp(r, t) and 
hence a local gauge parameter. The form of Aj reflects 
the north/south pole singularity in the spherical coor- 
dinate system, where ip is not well defined. Near the 
poles, we may choose a gauge in which 7 = =p(p, which 
renders the gauge fields well behaved either at the north 
or south poles {6 — or tt), respectively, but not both. 
It is thus necessary to use different gauges locally in re- 
gions where the texture passes through both north and 
south poles. We emphasize that such singularities have 
a purely mathematical origin arising from our choice of 
coordinate system, and may occur where the texture is 
perfectly smooth. 



The product of the transverse components in Eq. ( 14 ) 
is a gauge-invariant second-rank tensor: 



= sin2 Odiipdjip + diOdjO 

+ i sin 9{di9d.j(p — djipdjO) 
^Gij + iJ^ij ■ 



(15) 



The real part Gjj is a kind of metric in spin space, which 
will not appear in the final results of this paper.ES By 
gauge invariance, only the Berry curvature, i.e., the curl 
of the Berry gauge field Aj = —Aj, appears in any phys- 
ical quantities, 

Tij = ImA+Aj = AJA^-Af A} = diAj-djAi . (16) 

In the rest of the paper, where necessary, we will de- 
note the Berry gauge fields by Af — — (cosSp =F l)9/(/3p 
when well defined on the north/south pole, respectively. 



Geometrically, the Berry curvature gives the solid an- 
gle spanned by the spin texture mp(r,t) per area 
in the {x^,x'-') plane. Nonvanishing Berry curvature 
means that particles acquire a phase-space Berry phase 
{s/2)fdxiAi = sQ./2 over a closed trajectory, which 
modifies their transport. All the phenomena we will in- 
vestigate in this paper may be traced back to this phase. 

Equation ( 11 ) may be viewed as an expansion of Eq. ([3]) 
in powers of h. However, the separation into its classical 
0{hP) and quantum 0{fi) part is not manifest because 
of its matrix structure, which represents the dynamics 
of quantum-mechanical internal degrees of freedom. In 
the next section, we will derive a kinetic equation for the 
distribution function projected on each band defined by 
Eq. (|8|, which systematically captures all 0{h) quantum 
corrections to the classical Boltzmann equation. 



B. Decoupling 



The diagonal part of kinetic equation (ITT|) reads 



d'epsd, 



Fp,(a„A,;e„n„n) = 0, (17) 



where Fps are the 0{h), "anomalous" terms in the Boltz- 
mann equation. We denote transverse (i.e., x,y) vector 
components of the distribution (in the local frame) by 
n. In the formulae throughout this paper, the partial 
derivatives di acts only on the symbol to its immedi- 
ate right, so that in expressions of the form diAdjB, 
di acts on A only. At this point, the anomalous terms, 
which are shown explicitly in Eqs. (D2) and (D3| of ap- 
pendix [D] may not appear manifestly gauge invariant 
(while they certainly should be) ^ Let us now decouple 
the longitudinal and transverse components in an adia- 
batic approximation, which will result in a closed, gauge- 
invariant equation for the diagonal distribution functions. 
The decoupling procedure may be organized in the fol- 
lowing way. Suppose that we can solve the transverse 
components in terms of the longitudinal in a gradient 
expansion to the (p — l)th order in space-time deriva- 
tives: n = n(") -I- n(^) + . . . -I- n^P~'^\ By substituting 
n in -Fps, which is at least first order in gradients, and 



droppig any (p + l)th-order terms in Eq. (171, we arrive 
at a pth-order equation for Upg. We will carry out this 
procedure to second order (p — 2), consistent with our 
initial gradient expansion in Eq. (11). 

From the off-diagonal part of Eq. (11 1, in regions where 



A ^ 0, one may readily find the transverse components 
to first-order space-time gradients [see appendix [D] for 
details] : 



-^d.ep 



din. 



(18) 



where we have defined the phase-space particle-number 
and longitudinal spin densities, n — -\- and Uz = 
n-f — n^, respectively. The expression includes terms con- 
taining the expansion parameters of the adiabatic ap- 
proximation: {hvp/ A)dr and {h/A)dt, as well as a term 
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0(1/ A°) that shows that interband coherences may be 
generated by gradients of the distribution, irrespective 
of the size of the gap and even in the absence of exter- 
nal fields. We will see an example of this in the Dirac 
equation in Sec |III A] 



Let us further consider why, as shown in Eq. ( 18 ), inter- 



band coherences (encoded by the transverse distribution 
function) are proportional to the transverse gauge fields. 
Recall that the Berry gauge fields are the matrix elements 
Af^ = i{m.,s\di\m,s'), where the local eigenspinors of 
our spin bands are given by |mp(r,t),s) = Up{r,t)\z,s). 
The transverse gauge connection determines the overlap 
between spin up and down states of nearby eigenspinors, 

(m(a;, t),s\m{x + dx,t + St),s') = -iAf'Sx^ . 



Consider the terms 0{nz/A) in Eq. (18l, in the presence 
of an electric potential in the energy. The particle spins 
are locally misaligned with the texture m by an angle x 
given by^ 



tanx 



A' 



where the applied electric field is E = — c)re/2, and we 
defined group velocity Vp = dpe/2. {I runs over d spatial 
dimensions.) A finite x stems from the slight misalign- 
ment of electrons spin with the local quantization axis m 
due to (i) drift in momentum space with velocity E, (ii) 
drift in real space with velocity Vp, and (iii) dynamics 
of the texture. In particular, a locally spin-up electron 
will have amplitude on the nearby down band propor- 
tional to A'l . [See Eq. ( [ll| ) for a microscopic expression 
for A'l .] For a near-equilibrium distribution function 
which depends only on energy, ?T.ps(eps), the ©(l/A") 
term reads deUpsdieA^ . Evidently, in this case, it also 
stems from electron drift, but comes from a different 
part of the phase space. For a Fermi-Dirac distribution, 
this term originates from electrons on the Fermi surfaces, 
while the 0{nz/A) terms come from the momentum re- 
gion between the Fermi surfaces. 



Substituting Eq. ( 18 ) in the anomalous terms in 
Eq. (IT7|, we find 



1 



1 



djUs 



(19) 



Here, the Berry curvature terms are expressed as 



J^ij = z • A/ X Aj 



recalling that index I runs over d spatial dimensions. We 
have defined the fictitious charge Qs = —sh/2, s = ±1 for 
the t,i bands, respectively. Equation (191 constitutes 



a central result of this paper. It is clear that these are 
0{h) corrections to the Boltzmann equation (17). In the 



following sections, we will omit h for convenience (setting 



H = 1). In the rest of this paper, we will explicate Eq. ( 19 ) 
and apply it to specific examples. 



The first line in Eq. ( 19 ) represents corrections to the 
single-particle equation of motion. Including only these 



terms, the longitudinal transport equation (17) would 
read 

dtris - [{d^es + qsJ'ti)J'^ + qs^^esJ''^] d^Us + . . . = . 

(20) 

We first note that there is a correction to the single- 
particle energy, Cp^ — e^s + Se^s, where 



(5er 



ApJ- 



(21) 



which, in particular, is related to the magnetic moment 
of a semiclassical wave packet [see Eq. (l33l)].13SlThe other 



Berry-curvature corrections introduce Hall-like terms to 
the single-particle equation of motion, which we define so 
that the terms shown Eq. ( 20 ) constitute the phase-space 



advective derivative of the distribution function, 



Dins 



{dt + xid.j)ns 



We t hus id entify phase-space velocites 
(r.,p,),l™l 



^3 



is = <9pes 



Ps = -dr^s - QsiJ^tr + Trpidr^es ~ FrriOp^ts) ■ (22) 



As noted in Refs. fTTHTSl these equations of motion are 
noncanonical and thus appear to violate Louville's theo- 
rem. Indeed, the phase-space velocity has a finite phase- 
space divergence, 



dixl 



(23) 



From classical considerations of a two-component fluid, 
we would expect the projected kinetic equation to read: 



dtrir- 



V • {nps±s 



where we denote the phase- 



space gradient by V = (9r,9p), and . . . represents pos- 
sible interband terms. It may be seen, using the identity 
hd^Fij — djF (see appendix[E|), that terms proportional 
to to Eq. ( 23 ) indeed appear in the second part of the sec- 
ond line in Eq. ( |19[ ), however, multiplied by the longitudi- 
nal spin density which couples the bands. "^"^ To express 
our sKE as a kind of phase-space "continuity equation," 
we write — nz/4 = (qsUs + q-sn-s)/2 = qgUs — qs'n/2 in 
the expression for Fg, then our sKE reads 



dtHs + di{nsxl) = 



qs^{df + d'esd,) 



F, 
(24) 

where we have introduced notation for the relative band 
velocity: v\ = d^A. 

Both sides of the equation contain terms that violate 
Louville's theorem. Consider flrst the LHS. Since the 
phase-space divergence is in fact the advective derivative 
of qsF (up to the order of our gradient expansion) 



d,xl^qs{dt+d'esd,)F 

^ qs{dt + xld,)F + Oidl) 



QsDtF, 



(25) 
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we are lead to the following explanation for physical ori- 
gin of the term Ugdixl on the LHS. In a semiclassical 
description, spinful particles occupy wave-packet states 
labeled by momentum and position, as well as a spin- 
coherent state labeled by spherical angles on the unit 
sphere. A particle is thus specified by a set of continuous 
coordinates, (r, p, il), which we will call the semiclassical 
state space. However, in our adiabatic approximation, we 
retained only part of the (matrix) distribution function, 
nps(r, t), which denotes particle occupation per solid an- 
gle SV, subtended by the texture mp(r,t) in a volume 
Yli^riSpi. Since this solid angle is determined by the 
phase-space Berry curvature, being for example in one 
spatial dimension Sil{p,r,t) = J^rpSrdp, the changes in 
the distribution function due to the modulation of this 
spin space volume along the phase-space particle trajec- 
tories appear in the transport equation ( 24 ) as an advec- 
tive derivative of T i 



In Ref.[TTJ it was assumed that satisfies Eq. (24) 
with RHS equal to zero, in which case one could re- 
cover Louville's theorem for a rescaled distribution func- 
tion /s defined by = 2?s/s!^ where Vg satisfies 
d,x\ = -{l/Vs)DtVs, so that the LHS of Eq. ([24| reads 
VsDffs- Since Vg satisfies diil 



-DllnVs, from Eq. 



(251 we find 



V. 



1 - <ls^ , 



(26) 



Since we are neglecting terms that are cubic order in spa- 
tiotemporal derivatives at the level of the sKE, and, to 
the lowest order, Dt ~ 9^, we should keep only to 
linear order in gradients, as in the approximation above. 
This expression agrees with the one quoted in Ref.fT^ 
In Ref.[ITl Vs was included as part of the phase-space 
measure representing a modified phase-space density of 
states for band-projected electrons. Here, for our two 
band, continuum model, we have given an explanation 
for its origin. It accounts for the (pseudo)spin degrees 
of freedom in an approximation in which the quasipar- 
ticle spin dynamics is effectively constrained on a sub- 
manifold of the total state space. This submanifold, 
which we will call the projected state space, is a hyper- 
surface {r, p, il(r, p, i)} determined by the spin texture 
m[f2(r, p, i)]. The projected state space has a local, dy- 
namical volume proportional to J^(r, p, t) and is thus 
curved. We note that this curvature vanishes, J- = 0, un- 
less mp(r,<) has both momentum and real-space deriva- 
tives. 



Now consider the terms in the RHS of Eq. ( 24 ) , which 



contain interband terms. These terms are to be expected 
because the wave packets in the projected state space 
are not pointlike, but occupy a finite minimum spin solid 
angle All corresponding to a finite minimum volume 
Y\i Ar;Ap; ^ in phase space required by the uncer- 
tainty principle. However, due to the nonorthorgonality 
of the spin-coherent states, the spin up/down (along the 
texture) wave-packet states have nonzero overlap, result- 
ing in interband couplings in the transport equation. The 
overlap amplitudes are proportional to the area AJl ^ 



consistent with the fact that the interband terms are pro- 
portional to J- . 

It is clear that these effects are generic to projected 
kinetic equation in multiband systems. Since the decou- 
pling procedure is technically much more difficult in this 
case, let us sketch out qualitatively how our two-band 
result might be generalized. In a low energy, k • p ex- 
pansion about a point (in the Brillouin zone) of high 
symmetry, electronic quasiparticle transforms under an 
irreducible spinor representation of the crystallographic 
point group, which in some cases may be taken to be a 
higher-spin group, as in the case of the Luttinger Hamil- 
tonian. The multicomponent electronic wave function 
may be viewed as an amplitude on a group manifold (in 
the representation space of the symmetry group). There- 
fore, we may specify quasiparticles as in the two band 
case, except now il denotes coordinates on a higher di- 
mensional group manifold, and Eq. (Ill, which governs 
coupled orbital and internal dynamics, may be viewed as 
quasiparticle motion in (r, p, fi). 

For homogeneous spin-orbit couplings, the eigenfunc- 
tions "i/js = Ms(p)e'P''' of the k • p Hamiltonian H(p) 
defines a set of momentum-space spin textures. When 
the spin-orbit field is inhomogeneous, we can still de- 
fine a local, plane-wave-like (overcomplete) basis ips = 
Ms(r, p)e'P '', where Us(r, p) is an eigenvector of the 
Wigner-transformed Hamiltonian, and build wave pack- 
ets out of these states. Thus, the situation is quite similar 
to the two-band case, and we expect similar arguments 
pertaining to a curved state space and interband cou- 
plings due to the nonorthogonality of semiclassical co- 
herent states to hold. 



C. Hydrodynamics 

It is natural to expect a continuity equation for the par- 
ticle density, p — Tr[np] — X]ps "-ps^ which follows by 
taking X]ps '^^ '^^^ sKE's, resulting in an equation, 



dtp + dr-y 



= 0, 



from which we may extract the curl-free part of the par- 
ticle current. The anomalous, 0{h) contribution to dtp is 
given by the terms in Eq. (19), which, upon summation, 
may be written as [cf. Eq. ( |D13] )] 



d'^p 



ps \ Sis ^'^p-^ 



-a J" 



(27) 



where we have denoted the Berry-curvature corrections 



proportional to Qs in Eq. (22) by (5rs,<5ps) (these are 
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the corrections not including the shifts in band energies 
diSes). The momentum-space integration extends to a 
cutoff Pa above which our low-energy effective theory is 
no longer valid. 



From the first line in the second equality of Eq. (27) 



and including the "normal" part of the Boltzmann equa- 
tion Eq. (If 7p, we identify the particle current 



- QsiJ'tp + J'ppidri^s — J'pridp,€s)] . (28) 



Note that only part of the energy correction in Eq. (2f ) 



and the corresponding contribution to group velocity en- 
ters into this current. This resulted from cancelations 
which occurred in tracing the RHS of Eq. (24). 



The second line in Eq. ( 27 ) appears as an anomalous 



particle source coming from the boundary in momentum 
space (invoking the divergence theorem). Thus the ... in 
our continuity equation reads, 



^i72^^-K^ 

Qs {J' tr ^ J' cpi ^ri J' rri ^pi ^s) J" ; (^'^j 

where 1 is the normal to the bounding surface at p/^. 
This term represents particles flux coming from outside 
the momentum-space region where our semiclassical the- 
ory is valid. In principle, if we know the microscopic 
theory beyond the cutoff and could solve for the distri- 
bution function there, we may use it as an input in our 
low-energy theory, providing a momentum-space bound- 
ary condition for our distribution function that is found 
below the cutoff. In the following, for simplicity, we will 
take the cutoff to infinity. It is well known that iso- 
lated degenerate points (band crossings) where Ap = 
are monopole sources of Berry curvatur d^ l ^ l ^^ l (see ap- 
pendix |e]). These singularities occur because it is not 
possible to choose a spin frame [defined by the SU(2) ro- 
tation in Eq. ^] that is continuous, near topological de- 
fects in the spin textures.'^ This occurs, for example, at 
the origin of the hedgehog texture in the Hamiltonian of 
Eq. (49 1. Strictly speaking, our expression for the trans- 



A. 



verse component in Eq. ( 18 1 is not valid at these points 
(since A = 0).^SlQenerally, we will have to exclude these 
points by imposing a lower bound in our momentum inte- 
gration defined by an infinitesimal surface bounding the 
singularity. The boundary conditions for rips at these 
"inner" boundaries may be found by solving the exact 
quantum-mechanical problem near the crossing point. 

We note that since all terms oc A~^ canceled out in 
the final transport equation, we may take the dc limit: 
w — >■ 0, A — >■ (in the prescribed order). Thus in some 
instances, one may regulate monopole singularities by in- 
troducing a gap in the Hamiltonian, compute the current, 
and then take the gap to zero. However, this procedure 
will in general introduce some ambiguity in the final an- 
swer. For example, it may depend on the direction of 



the texture at the point where the gap is taken to zero. 
A well-known example of this problem is the case of the 
massless 2D Dirac fermion,!^ where the vacuum current 
depends on the sign of the mass used to regulate the di- 
vergences for the massless fermion. Below, we will work 
out the massive case in our semiclassical approach. 



D. Coupling to electromagnetic fields 

Consider a situation where the real-space texture in- 
homogeneity stems only from minimal coupling to vec- 
tor and scalar gauge potentials, a(r, t) and (/)(r, t), in the 
first-quantized microscopic Hamiltonian, 

n {-idr, r,t)^n {-idr - a(r, i), r, t) + <j>{r, t) . (30) 

The electron charge is absorbed here in the definition of 
the gauge potentials. 

From the Wigner transformation of this Hamiltonian, 
the semiclassical energy becomes 

•^p fp-a(r,t) + 0(i", t) = ek(r, t), 

where k(r, t) = p — a(r, t) is the kinetic momen- 
tum. The energy and the canonical momentum are 
not gauge invariant. To make the transport equation 
manifestly gauge invariant, we must express the dis- 
tribution function in terms of the kinetic momentum, 
= "'k+a(r,t)(r, i) = "■k(r,t). Note that 9k = 5p, but 
the space-time derivatives 9^ = {dr,dt) in the transport 
equation are taken with fixed canonical momentum p. 
We have to take account of the implicit space-time de- 
pendence of kinetic momentum by 



9^|p = 9, 



9pk • 9k = 



dnSL ■ 9k : 



(31) 



where it is implicit now that 9^ is taken with fixed k. 
We will also impose a gauge invariant momentum-space 
cutoff k\. The transformation to kinetic momentum may 
be viewed as a phase-space coordinate transformation, 
under which the gradients transform as 



dr^ps — 9reks — vi^sjdj-ai , 

9r71p 9r?lk - 9fe,nk9ra; , 



(32) 



Here, the group velocity is denoted by Vks = 9kes- 
The transformation of the Berry-curvature corrections 
are given in appendix |Fj We find that = b • Bk, so 
that the energy shift (5ek can be interpreted in terms of 
a magnetic moment Alk, 



5ek 



-Mk b 



(33) 



where A4k = AkBk/4. The factor of 1/4 is a result of 
our conventions, where we have factored out 1/2 from 
the usual definition of the Berry gauge fields in the lit- 
erature. We show in appendix |H] that this agrees with 
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the magnetic moment derived more directly from wave- 
packet analysis?^ The gauge-invariant transport equation 
now reads: 



4 



{dt + vks • dr + fks • 9k) (b • Bk) 

(vk^ • 9r + Vk2 X b • 9k) ?^k-s = , (34) 



where we have defined the relative velocity Vk^ = Vk^ — 
Vk^ and Lorentz force f, = e -I- Vks x b. e = —dr(j> — dtSi 
and h = X a are the ordi nary electromagnetic fields. 
The phase-space velocities are^E^ 

i"s Vks + 9k(5ek - Qsiks x Bk 

= ^ksVks + 9k<5ek -qs[exBk + (vks • Bk)b] , (35) 
kg fs + Sts - (7s(fks X Bk) X b 

= Zk,f. + 5fs-g,s(b-e)Bk, (36) 
where we define Zks = 1 + q^b • Bk and 
Sfs — —drSek + 9k(5ek x b . 



We note that Eq. (36) is not simply in Eq. (22 1 trans- 



formed to kinetic momentum (which would not even be 
gauge invariant ) , as it was necessary to transform the en- 
tire sKE. In the second equalities for , kg , we see that 
the group velocity and electromagnetic force are rescaled 
and shifted. The shift Siks simply comes from the mag- 
netic energy and the additional group velocity due to the 
magneto moment in Eq. (33). The total 0{h) correction 
to the group velocity, if we write the first two terms in 
the second equality for in Eq. ( 35 ) as Vks + iJvks , is 



<5vks = gs(b • iBk)vks + 9k(5ek 
= -(7s (b • Bk)(vk 



A. 



svk.) - ^9k(b • Bk) . (37) 

The third term in the first equality in expression for 
is the well-known anomalous Hall velocity.Sl Furthermore, 



the last terms in the second equality of Eqs. ( 35 ) and ( 36 ) 



give an additional phase-space velocity in the direction 
of the phase-space "magnetic fields." 

Up to the order of our gradient expansion, in which 
we need to keep terms in Ts (which multiplies drUs) to 
linear order in EM fields, and in k^ to quadratic order, 
and defining Vks = dki^s, e^s = -drius - Sja (recall 
that Eks — Eks + S^ks), the equations of motion may be 
expressed as 



= Vks - gsks X Bk 
ks = Gks + i-s X b , 



(38) 



in agreement with the form of wave-packet equations of 
motion.-^ We emphasize, however, that we are retaining 
terms beyond linear electromagnetic respons^^ which 
are not in the usual wave-packet equations of motion, 
for example in the b • e term for kg in Eq. ( 36 ) , and in- 
cludes forces such as —d^Sey^ coming from magnetic field 
gradients. 



III. EXAMPLES 



Massive 2D Dirac fermions 



As an application of our transport equation and a 
check on our formalism, we consider the electromagnetic 
response of massive 2D Dirac fermions, which is relevant, 
for example, to the gapped surface state of 3D topological 
insulators. The Hamiltonian is given by 



Hd = ^ ^a[-iOr - a{r, t)]a + mv^az , (39) 



where w is a constant with the units of velocity (which 
will be set to unity henceforth) and o- is a vector of Pauli 
matrices. The particle/hole symmetric dispersions are 
conveniently expressed in terms of the relativistic ener- 
gies Ek = Ak/2 — Vk^ + m^, so that eks = sEk with 
corresponding group velocities Vks = s'k./Ek- In k space, 
the texture is a vortex (meron) centered at the origin, 
V'k — arg(k), with an out-of-plane component given by 
cos^k — m/Ek- The vortex polarity is given by sgn(m). 
The spin texture has Berry curvature 



Sk = -9k cos 6'k X dkipi, -^z . 

The magnetic moment is thus A4]i — m/2Ef.. Both Berry 
curvature and magnetic moment are localized near the 
origin (see Fig.[l]) and have a direction (normal to the xy 
plane) depending on sgn(TO).'^ 




FIG. 1; In dimensionless variables, a plot of the Berry cur- 
vature, nr^Bk (solid curve), and magnetic moment, mMu 
(dashed curve), as functions of k/m. 



Consider now applying a static, homogeneous magnetic 
field b = 6z. From Eq. (37), the correction to the group 
velocity is (5vks = (fem/2i?^)k. The equations of motion 
Eq. (|35l) and Eq. (pel) then read 



rs = s + 



bm \ k 
2El) Ek 



k, = 6r, X z . 



(40) 
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It is evident that the usual cyclotron motion holds, but 
the cyclotron orbits are not exactly traversed in the op- 
posite sense for the two bands because of the 0{b) cor- 
rection. We also note that to reach this result, we had 
to keep the 0{b^) term in the Lorentz force, which arises 
due to the magnetic field dependence of the group veloc- 
ity. The second line in ( 34 1 vanishes on account of the 



spatially uniform static magnetic field and the fact that 
fks is transverse to k. The third line reads 



B 



[bk ■ dy + 6^k X z • c'kj'^k- 



(41) 



For a spatially homogenous distribution (which is valid 
for an infinite plane), the 0{b) term vanishes, and for a 
rotationally invariant distribution in k, the 0{b^) term 
vanishes. Thus, in this case, the transport equation is 
simply solved by a constant distribution function. 

To see some dynamics, consider adiabatically turning 
on the magnetic field, and let b{t) be a slowly increasing 
function. We choose a rotationally-invariant vector po- 
tential that produces a uniform magnetic field along the 
z axis, a = 6z x r/2 = b{—y, x)/2, with the accompanying 
circulating electric field e — — 9ta = —bz x r/2. First, 
let us check that our equations agree with the Stfeda 
formula,'231 which is a general relationship between the 
Hall conductivity in a gapped system and the change 
in particle density p when a magnetic field is turned 
on adiabatically. For the rest of this section, let us re- 
store h, the speed of light c, and the electron charge by 
e = eE, b = eB/c. In 2D, the Stfeda formula reads 
axy = ecdp/dBz. 

For a homogenous distribution function in linear re- 
sponse, the transport Eq. p4|) reduces to 



e ■ 

2^^ 



(42) 



where is the unperturbed, longitudinal spin distri- 
bution that is zeroth order in B. The resultant change in 
particle density Sp for an adiabatic change 6Bz is given 
by 



ec 



SBz 



(43) 



where 



.(0) 



Bu 



The extra particle density comes from particle fluxes at 
infinity due to the circulating electric field that produces 
anomalous Hall velocity ths = — eg^E x Bu- This results 
in particles on the upper (lower) band to enter (leave) 
the system. See Fig. [2] One finds the corresponding Hall 
current 



^^r^srik, 

ks 



(0) 



u—E, X z . 

h 



(44) 




FIG. 2: The anomalous current (Gl I comes from particles in 



the hole band flowing out of the system due to the anomalous 
Hall velocity ((ill. 



Comparing with Eq. (43), the conductivity axy = ve /h 



satisfies the Stfeda formula. Note that this is a general 
dynamical result, and, in particular, we do not need to 
assume a Fermi-Dirac function for the unperturbed dis- 
tribution. If the unperturbed distribution function de- 
scribes the ground state given by n|^°J — Q^ep — Ckt) ^ 
6(eF - Cki), then 



An 



(45) 



where $s = / (PkQ{ep — e^s) B^ are the Berry curvature 
fluxes over occupied regions on each band.^^ As another 
check, we show in appendix[G]that v indeed gives the cor- 
rect ground-state current of 2D Dirac massive fermions 
computed from field-theoretical methods. 

Let us now consider our results in light of the well- 
known physics of the Dirac equation. In this case, 
the adiabatic approximation corresponds to the semi- 
classical limit, valid in the presence of weak external 
fields in which particle/antiparticle pair creation may be 
neglected.'^ Historically, the single-particle interpreta- 
tion of the Dirac equation showed some puzzling features. 
The velocity operator given by c& in 2D has discrete 
eigenvalues ±c [restoring units here for the elementary 
electron, for which v = c in Eq. (39)], which may seem to 



contradict the experimental fact that the electron's ac- 
tual velocity is much less than the speed of light. Futher- 
more, it has noncommuting components and, therefore, 
does not lend itself to a classical interpretation. Dirac 
first resolved this apparent paradox by demonstrating 
that the electron trajectory exhibits a trembling mo- 
tion called the "Zitterbewegung:"'^ The electron posi- 
tion moves with mean (group) velocity pc^ /Ep and fast 
oscillations with frequency of the order of the mass gap 
A — 2mc? /h and amplitude of the order of the Comp- 
ton wavelength Ac = h/mc. However, when the electron's 
rest-mass frequency and Compton wavelength are beyond 
the finite time and spatial resolutions of real experiments, 
only the mean velocity and position are measured. Our 
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transport equation correctly captures the transport cur- 
rent due to the mean velocity. 

By considering the current of an electron wave packet, 
one can show that the group velocity pc^/Ep conies 
from the positive-energy band, while the Zitterbewegung 
comes from mixing with the negative energy bandl^^l Fur- 
thermore, it can be shown that the minimum size of a 
wave packet constructed from only positive-energy states 
is (in the nonrelativistic limit) of the order of In 
other words, it is necessary to have a mixture of both 
positive- and negative-energy states to localize a wave 
packet in a spatial region smaller than Therefore, 
a semiclassical theory where the electron velocities are 
given by pc^/Ep requires a one-band description, which 
is possible only if the size of electronic wave packets is 
much larger than Ac, even in the absence of external 
fields. 

This requirement shows up in our formalism as follows. 
In the expression for the transverse distribution function 



n in Eq. (18), the term 



hdriuAp^ must be small for 



our approximations to hold. Consider the magnitude of 
transverse gauge field in a relativistic expansion in the 
particle speed, /Sp = |vps|/c — pc/Ep = sin0, 



n\Ap, I = h^smH{dp^)^ + {dp9y 



hcos e^ian^ d{dpipY + cos^ e{dp tan d f 
^cos0Vl + cos2 = ^^(l-/32)(2-/32) 
V2 



A, 



27r 



OiP'p) 



Xr 



(46) 



where we used tan^ — p/mc, {dpipY — Thus we 

require |fip/«| ~ h\drinApJn\ ^ Xc/l ^ 1, where I is 
the lengthscale for spatial gradients of Up. That is, the 
distribution function must be smooth on the scale of the 
Compton wavelength, consistent with the size require- 
ment on nonrelativistic wave packets. Note that in the 



ultrarelativistic limit, the right-hand side of Eq. (46) ap- 
proaches the de Broglie wavelength h/p, which then has 
to be much smaller than I. This is always required in 
the gradient expansion, even in the absence of interband 
terms. 

The microscopic currents associated with the finite size 
of the electronic wave packet lead to an appearance of the 

EkBk/2. In fact. 



magnetic moment [see Eq. ( 33 )] A4k 



the anomalous velocity may be understood from the rel- 
ativistic physics of this magnetic moment, as foUows.l^ 
Under an applied E field, in the frame of an electron 
moving with velocity Vks = shc^/Ek, by the Lorentz 
transformation, there is a magnetic field given (for small 
velocities) by B(vks) = — Vks x E/c. The magnetic en- 
ergy in the electron rest frame then reads 



he 



he 



(5eks = ^Mk-VksXE=ysk-ExBk, (47) 

which in the lab frame may be interpreted as the energy 
of an effective electric dipole Pks = (^e/c^)vks x Alk- 



The contribution to the group velocity d^Se^s due to this 
extra energy is exactly the anomalous velocity. 

Finally, we point out some surprising features in our 
band-diagonal transport equation that call for further in- 
vestigation. Equation (34) (where the second line is pro- 
portional to riz) shows that the Stfeda formula does not 
hold separately on each band because half of the expected 
intraband flux goes into interband fiux, which seems to 
suggest the following phenomenon. Suppose the mag- 
netic field smoothly falls off to zero in an outside region 
far away from the origin, where there are well-defined 
energy eigenstates. Due to the anomalous Hall velocity, 
particles on the upper (lower) band flow in (out) of the 
sample. Inside the region with magnetic field, particles 
are transferred between bands. If we started with the 
lower band occupied in the outside region, and slowly 
turn on the magnetic field, we expect a pumping effect 
whereby some of the particles traveling into the origin 
would come out on the upper band (as defined by our 
WT basis). However, we note that the interband terms 
in the second and third lines of Eq. ( 34 ) that are induced 



by spatial inhomegeneities would also be present in this 
situation and need to be taken into account. The ques- 
tions of whether the adiabatic interband fluxes in the 
presence of a finite gap can be physically manifested or if 
one could contrive a transformation that eliminates such 



interband terms in our transport equation ( 24 ) altogether 
remain open. 

B. Faraday's law in momentum space 

As another instructive example, consider a slowly- 
varying gap m — m(t) of 2D Dirac fermions, which, for 
example, would occur due to a Zeeman splitting in the 
z direction in topological insulators. For simplicity, we 
assume here the magnetic field is static. Then the trans- 
port equation is 



{dt + rs ■ dr + iis ■ 9k)nks 



(48) 



where B — m/ E^. The time-dependent Berry flux causes 
an anomalous velocity [cf. Eq. (22)] 5vs — QsJ'kt = Qs^ks, 



where we have defined the momentum-space fictitious 
electric field which in this example is given by 



5k = dt cos 9dii(p 



m 



z X k . 



This anomalous velocity is transverse to the cyclotron 
orbits and gives a radial particle flow in opposite di- 
rections on the two bands, similar to the previous ex- 
ample, except that it is peaked at fc = m/\/2. This 
radial flow causes changes in particle density, deter- 
mined by the time dependence of the Berry curvature: 
B = TO(fc2 - 2m^)/El See Fig.jl It may readily be 
verified the momentum-space, fictitious electromagnetic 
fields satisfies Faraday's law, 9t6k -f 9k x £k = [See 



Eq. (E4) 
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FIG. 3: A plot in dimensionless variables of the fictitious elec- 
tric field, (m^/m)i£k (solid curve), and Berry fiux variation, 
(m^/m)Z3k (dashed curve), versus k/m. 



C. 3D Weyl fermions 

As another illuniimating example, we consider 3D 
Weyl fermions, which represent 3D band crossing points. 
In the presence of electromagnetic fields along the z axis, 
the Hamiltonian for a right-handed (RH) Weyl fermion 



n 



w 



vn ■ \—id-r — a(r)] 



(49) 



Here, a — a^^x -I- a^y is the vector potential associated 
with a magnetic field 6 > {Bz < for electrons) along 
the z axis, using the same gauge as in Sec. |III A] The 
electric field is given by = —dtaz(t). 



The Hamiltonian in Eq. ( 49 1 has semiclassical disper- 
sion Ek = and the Berry curvature of a monopole 
source at the origin: Bk = k/fc'^, which causes a 
momentum-space particle flux represented by th e b • e 



term for in Eq. (36). As discussed in Sec 



HC 



to ex- 
clude the singularity at the origin, we must introduce 
an inner boundary in the momentum-space flux integral 
in Eq. (29). The particle fluxes through this boundary 



produce an anomalous source of particle density, 



dtp = 



47r 



dSk • Bk Ukz 



S2 



BE. 



(50) 



Up to quadratic electromagnetic response, the integral 
in Eq. (50l depends on the zeroth-order (in E and B) 



distribution function n^^ on the boundary surface, which 
we take to be 5^. Assuming a ground-state distribution 
given by the Fermi level at zero energy, so that only the 
lower band is filled, n^^ = —&{eF — ^h), the integral 
gives 47r.'^ In units where h ^ 1, this results may be 
expressed as dtp = — (e^/47r^) B • E, which is known in 
particle physics literature as the (3 -I- 1)D Adler-BcU- 
Jackiw anomaly,^^ and, strictly speaking, cannot occur 
alone in nature because it violates charge conservation. 



Let us examine the source of this particle flux in the 



quantum mechanical solution of Eq. ( 49 1 . In the presence 
of a magnetic field along z, the motion in the x — y plane 
is quantized into Landau levels. By translational invari- 
ance, the momentum along the magnetic field remains 
a good quantum number. The spinor eigenfunctions are 
thus of the form ^'„(r) = e'^P'^ipp^n{x,y)- Since 

-Hiv*n(r) = e'^-'^nDim ^ Pz]iJp,n{x,y) , 

ipp^nix,y) are the eigenfunctions of the Dirac Hamilto- 
nian ( 39 ) with a mass Pz ■ The square of Hamiltonian 



(49) determines the eigenvalues (up to a sign), which are 
the relativistic Landau levels,!^ 



Ena = (TV^/2\eB\ 



Eo = vp^sgn(eS) . 



Here n — 1,2,... is a Landau-level index and a — ±\ 
labels the particle/hole branches. The Landau levels 
spread out into bands in the z direction, all of which 
are gapped, except for the "zero mode" E'q, which has a 
linear, chiral dispersion along the direction of magnetic 
field. Starting with an equilibrium distribution with the 
Fermi level at zero energy, if we apply an electric field 
along the z axis by adiabatically turning on az , particles 
are transported above the Fermi level by spectral flow, 
locally populating states near the momentum-space ori- 
gin. These particles originate at Pz = — oo (outside the 
momentum-space cutoff ).'35i 

The particle production rate may be computed from 
the rate at which the Fermi surface of the chiral branch 
changes, ppz = &Ez, taking into account the zeroth 
Landau-level degeneracy. The degeneracy of the Landau 
levels are given by ^/^q, where $ = |i?z^| is the mag- 
netic flux {A is the transverse area) and, in units where 
h — \ and c = 1, the flux quantum is $0 — 27r/e. Thus 
the density of states per area in the transverse direction 
is ei3/27r, and per length along z is l/27r. Therefore^^- 



dtp 



eEz e\Bz 
27r 27r 



fEzBz 
47r2 



(51) 



We may thus match this quantum-mechanical solution 
for dtp with the semiclassical result given in Eq. (50). 



We stress that our semiclassical analysis here involves 
response quadratic in electromagnetic fields and is, there- 
fore, beyond the scope of the wave-packet analysis found 
in the literature."^ 



IV. CONCLUSION 

In summary, we have derived a band-diagonal, semi- 
classical kinetic equation (sKE) for electrons in a two- 
band, (pseudo)spin-orbit coupled system, which takes the 
form of a Boltzmann equation [see Eqs. ( [l7| ) and ( 19 )] 
for a coUisionless plasma. In addition to the corrections 
to quasiparticle equations of motion, we find terms pro- 
portional to the distribution function that we attribute 
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to single-particle motion constrained on a curved state 
space. We also find interband couplings that represent 
coherences due to the nonorthogonal nature of the pro- 
jected Hilbert spaces. As a check on our formalism, we 
find our kinetic equation reproduces the well-known elec- 
tromagnetic response of 2D Dirac fermions and 3D Wcyl 
fermions. For Hamiltonians with less symmetry, more 
complicated inhomogeneous spin-orbit couplings, exact 
quantum-mechanical solutions are hard to find, while our 
sKE remains valid and a useful analytical tool. 

At first sight, the interband terms in our sKE seem 
puzzling, given that there are well-defined procedures 
for decoupling multiband Hamiltonian in the presence of 
weak external fields.** ^ However, these methods employ 
complicated quantum-mechanical transformations on the 
Hilbert space which generally do not constitute a simple 
change of basis. BlouniP used a mixed representation 
where operators are specified by both Bloch (or canoni- 
cal) momentum and position. By successive transforma- 
tions, that author derives an effective Hamiltonian for 
Bloch (Dirac) electrons in a magnetic field as a function 
of kinetic momentum i?(k), which is a power series in 
B. The transformations are functions of k and become 
nonunitary in the presence of a magnetic field because 
of the commutator: [ki,kj] = ieijkb^ . As pointed out by 
the author, qualitatively, these transformations amount 
to a "local" diagonalization of the Hamiltonian in phase 
space, similar in spirit to our approach. 

For the (3-|-l)D Dirac equation, one may apply the 
Foldy-Wouthuysen (FW) transformation,!^ to decouple 
the positive/negative energy bands in the presence of 
weak external fields in a power series of inverse mass 
1/m, resulting in the nonrelativistic Pauli Hamiltonian. 
The FW representation has the advantage that the op- 
erators retain their classical meaning. In the absence of 
external fields, one can show that the position operator 
r satisfies r = ip/E'p (for the positive/negative energy 
bands) , and the orbital and spin angular- momentum op- 
erators are separately conserved. However, these opera- 
tors correspond to the original Dirac operators evaluated 
at the mean position of the electron, which, as discussed 
in Sec. HI A is spread out in a region of (at least) the 
size of the Compton wavelength Ac. In the presence of 
external fields, in principle, the meaning of the operators 
also changes and it may be seen explicitly that part of 
the FW transformation becomes nonunitary. 

The key features of these methods are consistent with 
our formalism. Exact diagonalization of spin-orbit cou- 
pled Hamiltonian ■H(p) is possible only in the absence 
of inhomogenieties. When the bands are gapped, one 
intuitively expects that in the presence of smooth space- 
time inhomogeneities, quasiparticle motion could be con- 
fined to separate subspaces, which may be called field- 
modified energy bands .1^ One may indeed find represen- 
tations in which the Hamiltonian is approximately block- 
diagonalized in the field-modified band space, but at the 
expense of altering the Hilbert space with nonunitary 
transformation, which may render the subspaces of the 



field-modified bands nonorthogonal and change the phys- 
ical meaning of operators. This caveat is consistent with 
our explanation of the interband terms in Sec. |II B| In 
cases where the Hilbert space may be specified by contin- 
uous coordinates, we find it useful to visualize the field- 
modified bands as submanifolds with possible nontrivial 
interband orthorgonality relations. 

The Wigner distribution function gives us an exact 
phase-space representation of nonequilibrium quantum 
dynamics, where we can treat r, p as real numbers from 
the get-go, allowing us to diagonalize the semiclassical 
Hamiltonian ep(r,i) in a straightforward manner. Find- 
ing the exact correspondance between our matrix trans- 
formation of the distribution function and the aforemen- 
tioned transformation on the Hilbert space (where r, p 
are operators) seems to be a nontrivial task, which we 
will relegate to future work. Our approach based on the 
density-matrix formalism is quite general, allowing for 
the treatment of multicomponent fermions or bosons, and 
in principle allows for a straightforward incorporation of 
electron interaction effects in the spirit of Fermi-liquid 
theory. Lastly, our projection process described below 
Eq. (17) is valid up to higher orders, and, in particular. 



with some additional labor, one may carry out the gra- 
dient exp ansion to 3rd order, which may warrant further 
study.HSl 
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Appendix A: Gradient expansion 

The gradient expansion expresses the WT of a convo- 
lution of two kernels in terms of gradients of the WT of 
the individual kernels through the following formula.!^ 
Let a{h) be the WT of A{B), then WT of [A ® B] is. 



WT[A ®B]^a 



exp —{dr ■ dp - dp - dr 



ah + -d.ad'b - —dAad'd^b - 
2 8 



(Al) 



where . . . indicate higher order terms. For fermions, the 
gradient expansion is valid when the length scale ^ of 
spatial inhomogeneities is much longer than the Fermi 
wavelength A^^ '--^ I/q^f, where qp is the Fermi wave vec- 
tor. One may consider it as an expansion in q/qp, and 
q ~ 1/^ is the characteristic wave vector of the spatial 
inhomogeneities . 

In fact, in the gradient expansion applied to A® B — 
B ® A m Eq. aU even powers vanish if A and B 
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commute because pairs of contracted indices are sym- 
metric. Therefore, in the scalar (single band) case, the 
leading quantum mechanical corrections are third order 
and 0{1t?)I^ However, when Eq. ([s]) has nontrivial ma- 
trix structure, the second order expansion is required for 
capturing all semiclassical 0{h) terms. 



Appendix B: Covariance of the kinetic equation 



The transport equation (111 is covariant in the sense 



that it remains the same form in an arbitrary spin 
frame. Specifically, under a local SU(2) transforma- 
tion Up — > U^hpUp, Eq. represented symbolically 

as dfi/dt transforms as dfi/dt — ^ Wdh/dtU. The covari- 
ant derivatives transform similarly, with the appropriate 
change in the gauge potentials given by 



A, ^ Aid + iU^ d.,U . 



Expressing the matrix gauge potentials as U — 
exp(z77°f the infinitesimal transformation with 77° ^ 1 
is given by 



n n 77 X n, 
Ai Ai - ?7 X + a^T? . 



(B2) 



Evidently, the spin components of the distribution func- 
tion transforms in the adjoint representation of SU(2). 



equation. It will be convenient to first express the lon- 
gitudinal equation in terms of n and Uz- The trace and 



the fz component of Eq. ( 11 ) reads 
1 



where 



dtn - 
dtn-z - 

F = - 



{died'n + d,Ad'nz) + F 



0, 
0, 



(Dl) 



9,A(A* X z • n) -I- A(z • A, X d'-h) 



=z • At X n 

h 



AAi • hA 
1 



(9ie(z • A* X n) 



8 



z • (L>2 e X D^'J n) 



(D2) 



(D3) 



The vector part of the second covariant derivative, de- 



(•gj^^ fined by DiDjfi = didjn/2 + Df^n ■ t, reads 

D^jii —didjU + diAj X n + Aj X diU 
-I- Ai X djU + Ai X (Aj X n) . 



(D4) 



Let us first simplify the last term of Eq. (D3). The an- 



tisymmetric part of of the second covariant derivative 
Eq. (D4) vanishes, as one would expect since the SU(2) 



field strength is zero. Indeed, the antisymmetric part of 



the last term in ( D4 1 reads 



(Ar, X A -1 X n) 



-[(A, X A,) X n]'' 



Appendix C: Integrability condtion 

By their definition, the matrix Berry gauge potentials 
are pure gauge and so that their SU(2) field strength 
vanishes, which provides an integrability condition: 



5.A, 



Ai X A.I 



(Cl) 



This is a partial differential equation which relates the 
three component of the gauge field. Mathematically, in 
a simply connected region where Eq. (Cll holds, one 



can always find a frame in which the gauge potential 
vanishes.'SS Contracting indices in Eq. ( Cl ) gives an iden- 
tity, 



2d'A, 



-A' X A, 



(C2) 



which was used repeatedly to arrive at the final gauge- 
invariant expression in Eq. (D13|. 



Appendix D: Decoupling procedure 



where we define anti-symmetrization by A^^Bj^ = 
{AiBj — BiAj)/2. Therefore, the commutator of covari- 
ant derivatives vanishes 

[Di, Dj]n — [diAj — djAi + x Aj] x n = Fij x n = . 

(D5) 



where we've used the identity Eq. (Cll. The remain 



ing symmetric part of the covariant derivative Eq. ( D4 ) 
reads, 

D(^^Dj)n — df^n + ^(.jAj) x n + 2A(,j x dj^n 

+ A(, • n)Aj) - (A, • A^> , (D6) 

where we defined the symmetrization symbol A(^iBj) = 
{AiBj + BiAj)/2. Therefore, we will only need to com- 
pute that quantity 



(D7) 



of which we will only need the terms containing , since 
the terms containing n will be third order in gradients. 
Taking n — > rt^z, noting that in the rotated frame we 
already have e = Az, and that the vector product (D7) 



In this appendix, we show the computations that lead above involves only the transverse components of the see- 



to the anomalous terms in Eq. (191 in projected kinetic ond covariant derivative (D6I, we find 
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Pi^) = • { [Aa(,A,) + 29(,AAj)] X z - A^(,Aj)} 



{ 



Aa,n^)z • Aj X {d'k^ + a^'A* - A'Gj') 



1 



1 



(5,An, - Aa,n,)aJ J-/ 



(D8) 



In the third hne, we use identity (C2) and (E2|. This 



expression is gauge invariant, as it should be. Next, we 
want to solve for ii to first order gradients. Consider the 
transverse part of the transport Eq. (Ill to the first-order 
gradient expansion,!^ 



A 



1 r 



Dtn — — z X n 



0, (D9) 



H " 2 

where the transverse covariant derivatives are defined by 
Din = djii + AjZ x n + n^A/ x z , 



Die ~ AAi X z . 



(DIO) 



In regions where there is a finite gap, Ap ^ 0, Eq. ( D9 ) 
shows that n'-"-' = 0, so that the lowest nonvanishing 
order is n^^^. Therefore, to first order in Eq. (D9) we can 
drop terms ~ djii, Ajh, then 



A 



z X n'^^ 



1 r 

2 



n^Ai X z 
dke{nzA.^ x z) + d'^nAAk x z 



(Dll) 



Solving for n results in the solution Eq. ( 18 1 



In approximating Eq. (D9) with Eq. (Dll), we have 
neglected terms of 0{dtn), 0{hdp€ ■ A^,ndre ■ Ap), and 
0{dpe • driijdre ■ dpii), in comparison with the commu- 
tator term which is 0(An). Considering the near equi- 
librium case when n is localized on the Fermi surface, 



our approximation implies the limit [see also Sec. Ill A 



for additional constraint related to the second term in 



Eq. (18)] 



huj^ A, C > 



Ai^l^rel < A, 



(D12) 



where is the Fermi velocity, uj (^) are the charac- 
teristic frequency (length scale) of the dynamics (inho- 
mogeneities) of the system. The last condition on en- 
ergy gradients is a requirement on the size of driving 
electromagnetic fields. The limits above define the adi- 
abatic approximation. Also, same conditions as above 
with A — J- Ep {Xp ^ C etc.), for intraband adiabaticity. 



Substituting the expression for n in (Dll) into the 



anomalous terms of Eqs. (D2) and (D3), wc find 



+ \{n,d,A - Ad,n,)d'Fi' - \^Td,n,d'A . 



(D13) 



The anomalous terms in Eq. ( 19 ) are given by Fs — (F - 
sF,)/2. 



Appendix E: Bianchi identity 

The Bianchi identity for the Berry curvature reads 

diFjK + djFKi + OkFij = 0, I ^ J ^ K, (El) 

and follows directly from the definition of Berry cur- 
vature in terms of gauge fields if they are nonsingular, 
since d^jJ^jn] = d^jdjAx] = 0. Let I J = ij be phase- 



space indices and contracting identity (El) by multiply- 
ing Ejj I then 



d,.F' 



OkF = . 



(E2) 



We have used_(|E2|) repeatedly, for example, we trans- 

[dt+d'eA]F, (E3) 



formed Eq. ^ by, 

d,Ft' + d^r^dkes 



We have used this identity in the second line of ( 19 ) 



In the following, it will be necessary to introduce dif- 
ferential forms in order to use Stokes theorem in a space 
with more than three dimensions.'^ The Berry gauge field 
(connection) is a 1-form A = Ajdx^ and the Berry cur- 
vature is a two-form F = ^TjKdx'' A dx^ , where sum- 
mation over repeated indices are implied, and the com- 



ponents of these forms are given in Eq. ( 15 ) and the text 



following Eq. (16 1. Denoting the exterior derivative by 
d, we may write F — dA — d(cos 9) A d^p, which is the 
area 2-form representi ng s urface elements on the sphere. 
The Bianchi identity (El) states that dF = and rep- 



resents a set of homogenous Maxwell equations for the 
phase-space, fictitious electromagnetic fields defined by 



•Fqt — : 



where q G {r, p} is a 3D vector in phase space which may 
have indices in both r and p. In differential forms no- 
tation, Eq. ( |E1[ ) are the components of a 3-form, which 
may be expressed in terms of the fictitious electromag- 
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netic fields, 



1 



df -diFjK dx' A dx-' A dx"^ 

qe{r,p} 



- {dtBq + X fq); A dq"" A dq'' } ^ 0. (E4) 



Here (Imn) runs through three dimensions. The first 
term represent the absence of monopoles and the second 
is the phase-space Faraday's law for the fictitious electro- 
magnetic fields. 



However, near crossing points, (El) needs to be recon- 



sidered, since point degeneracies of the Hamiltonian are 



monopole sources of Berry flux^ and Eq. (El) assumes 



the absence of such singularities. Consider first the Berry 
flux over a sphere containing the crossing point in phase 
space. Modulo 47r, it is an integer topological invariant 
known as the Chern number: <i> = F = 47riV, N e 
''^2{,S'^) — Z, and represents the winding number of the 
mapping between spherical surfaces from (r, p, i) space 
to spin space defined by the texture, mp(r, t) : S^, 
analogous to the skyrmion number in ferromagnets. Us- 
ing Stokes theorem in the form df ~ Jg^^^ F, with 
dM = being the surface of the sphere M, we have 



$ / df ^ Attn . 
hi 



(E5) 



If we take the limit of M being infinitesimally small and 
containing a topological defect of the texture at q^; in a 
3D q space, only the first term of df in Eq. (E4| con- 
tributes to Eq. ( E5 ) , which implies thalP 



(E6) 



Here we have used the fact that near a crossing point q^, 
one may make a linear expansion of the gap vector 



[(q - q<i) 



9q]Aq 



Y^{q' - q^^)e,, 



thus the texture has winding number = ±1, with the 
sign determined by det(e"), where e° = 5iA°(qd) are a 
set of basis vectors defined by the gap gradients. The 
Weyl equation in Sec. |IIIC| is an example of a monopole 
texture. 

If, due to time-dependent parameters in the Hamilto- 
nian, the monopole position changes in time, q^ = qd(^), 
the current associated with the monopole motion will 
appear as a source term in Faraday's law. To see this, 
consider the plane defined by = (;^(io) in the 3D q 
space. Near t^, we may approximate 

q' - q\{t) ^ {q' - q'^{ta)) + q'a{to){t - to), 

Substituting this expression in the expansion above for 
the gap vector, it is evident that there is a monopole 
texture located at {q^,q2,to) in {q"^,q^,t) space, and 



that the basis vector in the time direction is given by 
et = q^ei. The Berry fiux over an infinitesimal surface 
containing t his monopole in {q"^,q",t) space, computed 
using Eq. ( E5 ) , will have contributions only from the 
second term of df in Eq. (E4), thuJl^ 



Imn 



5{t ~ to)S{q^^ ~ q7{to)W - qlito)) ■ (E7) 

The sign here is determined by sgn[(7^ det(e°)], and LHS 
of Eq.(E7l is evaluated at = (7^(<o)- In terms of 
monopole motion, this nonzero winding number of the 
texture comes from the fact that, on the — Qdito) 
plane, the texture A(q™, q", (7^(to)) is a vortex whose 
polarity, given by sgn(q^(io) — Qd{t)), switches sign dur- 
ing the course of time as the monopole passes the plane 
a.t t — tQ. To express the RHS of Eq. (E7) in terms of 
the monopole current, we note that the linear expansion 
qli{t) near to may be inverted to yield a mapping i((7^), 
which we may use to change variables: d{t{qlj) — to) — 
Wdl^ildi^) ~ 9d(*o))) and that on account of the delta 
function 5{t — to), we may evaluated the LHS of Eq.(E7) 
at = Then, 



-9q X £q](io) = sgn[det(en]4^qrf,53(q- 



q™(io))- 
(E8) 



Si milarly , we note that the 2D Dirac equation (See 
Sec HI B ) with a time-dependent mass m{t) oc t would 
give a source term in Faraday's law. 

Lastly, the magnetic flux through an open 2D hyper- 
surface S may be converted to a contour integral of the 
vector potential around the boundary dM, provided that 
we use the appropriate gauge. 



i dl'Af 
Jas 



[cos9{l)Ti]di(p{l)dl^ . 

(E9) 



The integral depends only on the path [9{l),ip{l)] traced 
out by the texture in spin space space (which can 
be many closed loops) as the phase-space contour dM 
parametrized by I G [0, 1] is traversed. Except when the 
cone angle 9 — 7r/2, this flux is generally not quantized. 



Appendix F: Berry curvature in general coordinates 

Under a general phase-space coordinate transforma- 
tion x^ — 2;*(x), the energy and the distribution func- 
tion transforms as a scalar e{x) = e{x{x)) = e{x), while 
derivatives transforms as a covariant phase-space vector 



y A 

dx^ ^ dx' 



- dxi 



(Fl) 



Since the Berry curvature is a product of derivatives, it 
transforms as a second rank covariant tensor 



= (A^ J-A), 
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where the bar denotes the Berry curvature with deriva- 
tives w.r.t. the coordinates x. In particular, when trans- 
forming to the kinetic momentum for electromagnetic 
perturbations, the spatial derivatives transforms as (31), 



resulting in the following transformations of the Berry 
curvature 



= diSi ■ dta X Bk 



[dja X Bk) 



^P^P, = -^fc.fej = ejji(Bk)( Jv.pj = {dia X Bk)j , (F2) 
while T transforms as 



E 



_9a 

Or, 



X Bi 



daP_ 



B 



b-Bk 



(F3) 



The simplest example of the transformations in Eq. ( F2 ) 



is in the presence of an electric field in the vector po- 
tential gauge. Then £p = Tpt = e x Kk represents 
the anomalous Hall velocity. This may be contrasted 
with the scalar potential gauge, in which the kinetic and 
canonical momentum are the same. 

In the presence of spatiotemporal inhomogeneities in 
the spin textures other than gauge potentials, the Berry 
curvatures will have space-time dependence not related 
to kinetic momentum, and there will be additional terms 



in the equations of motion Eq. ( 35 ) and Eq. ( 36 ) . The ad- 



ditional terms come simply from the Berry curvatures in 
Eq. (22 1 with space-time indices, with p — > k. For exam- 



ple, if the texture is due to a ferromagnetic exchange field, 
the Berry curvatures J>r and J>f are fictitious electro- 
magnetic fields that are known to mediate spin-transfer 
torques in itinerant ferromagnets.*^ 



Appendix G: Parity anomaly 



For the example in Sec. |III A[ it is instructive to com- 
pute V when the Fermi level lies in the gap, which gives 
the well known half-integer conductivity of 2D Dirac 
fermions.'22l Although the Berry flux in Eq. (45 1 is con- 



centrated near the origin, up to a sign they depend 
only on the texture far away from the origin. To em- 
phasize the topological nature of the number v, we ap- 
ply Stokes theorem to convert integrals to Fermi sur- 
face Berry phases which determines the solid angle sub- 
tended by the texture, taking care to use the appro- 
priate gauge, ^k = —(cos 0k — sgn(m))9k<i5k. Then 
V = — (l/47r)/^k • d\i — — sgn(m)/2, in agreement 
with the formula well known in particle physics litera- 
ture for anomalous charge current associated w ith th e 
parity anomaly of the massive 2D Dirac equation,'2iEZl 



(j'^) =sgn(m) — e'^"^F„;3. 



(Gl) 



In the formula above, is the electromagnetic field 
strength tensor, F^^ = dy,A^ - d^^A^, = ($,-A), 



Foi = Ei, Fij = —eijkB^ , the expectation value is taken 
in the vacuum state, which in the single particle picture 
is the "Dirac sea," with all negative energy levels occu- 
pied and thus corresponds to Fermi level in the gap. One 
should restore units in ( Gl I by setting 27r = /i to com- 



pare with our result for v. The dependence on sgn(m) 
is a signature of the infrared singularity at k = when 
TO = 0, as discussed in Sec. |II C[ The 2D Dirac mass 
term explicitly breaks parity, defined by inversion of one 
spatial coordinate. In field theory, even when to = in 
the Hamiltonian, gauge-invariant regularization of the ul- 
traviolet divergence may violate parity (for example, in 
the Pauli-Villars regularization), leading to the parity- 
violating anomalous current .'SSl 

In field-theoretical methods, one calculates the anoma- 



lous current (Gil by the coupling of fermionic vacuum 



energy to an external vector potential, represented by 
a fermionic functional determinant. While the result 
is precise, its physical origin is somewhat mysterious. 
The semiclassical approach gives an intuitive but rigor- 
ous single-particle picture, in which the vacuum current 
simply represents the flow of quasiparticles. 



Appendix H: Comparison with wave-packet theory 

According to wave-packet analysis,'^ the correction to 
the single particle energy due to inhomogeneous pertur- 
bations is given by. 



5es = Im (driUslHcir^t) - es{Tc,t)\dp^Us) 



(HI) 



where i?c(r, t) is the local Hamiltonian, with the position 
r evaluated at the center of the wave packet, and £s{r,t) 
is the band s local eigenvalue. We emphasize that in 
the wavepacket approach, r is a c-number parameter of a 
quantum Hamiltonian, whose complete set of eigenstates 
is known. The local Hamiltonian is the zeroth order term 
in an expansion of a Hamiltonian with smooth inhomo- 
geneous perturbations. The perturbations to flrst order 
in gradients has the form 5H ~ (f — r) • drli and it's 
expectation value in a wave packet state gives the en- 
ergy correction in first order perturbation theory. From 
the semiclassical transport equation approach, the local 
Hamiltonian is equivilant to the semiclassical quasiparti- 
cle energy ([7|, 

Hc = ip{r,t) 

To compare with formulae from wave-packet analysis, we 
first define the local eigenkets by a rotation of the lab 
frame spin eigenkets pointing in the ±z; direction: 

\ups{r,t)) = Up{r,t)\i;s}. 

They satisfy the local eigenvalue equation: 

[i7,(r,t)-e,(r,t)]|up,(r,t)) = 0, (H2) 
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where is given by (10). For brevity, we will denote 



"lab" frame spinor basis by |s), s = ±1. The matrix 
elements of the gauge fields are given by 

[A, • tY^' = ir' = i{s\U^dMs') = i{us\d,u,,) , (H3) 
where we've omitted phase-space subscripts for brevity. 



From (HI), we find 



Se, = ^(s|a,C/^[Amp • f - sA/2]d'U\s) 

= -{s\d,U''U[Ar, ~ sA/2]C/^a*L/|s) 
2i 



1 

2i 



A ( \- 




- ^ 



(H4) 



where J- = Ti 12. Note that the energy shift is indepen- 
dent of s. 



Appendix I: Relation of gauge fields to velocity 
operator 



We note that the off-diagonal gauge field is related to 
the phase-space velocity operator, — diH. Differen- 



tiating the eigenvalue equation (H2) and projecting on 
\usi) for s s', we fi nd an equation for the off-diagonal 
gauge fields [cf. (H3)]: 



{us\diH\us 



(11) 



Consider going into the frame of the electron moving 
in phase space, then the phase-space dependence of the 
texture appears to the electron as dynamics of a local 
phase-space magnetic field. Therefore, to make heuristic 
comparison with the quantum adiabatic theoremj^^S! we 
take di — )• dt in Eq. ( [TT| ) . Then the theorem shows that 
the amplitude for transitions to other states is given by 
Af^ /A in agreement with Eq. (18). 
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In particular, AiA^ = 0. Note that since J" = -1, the 
inverse is J = — J. 
''^ Qij = ReA'I^A^ is related to the field strength by 



IJ ~ £IJ 



IJ ■ 



where we consider the metric tensor Qu as a matrix in the 
(xi,x,j) plane and defined its determinant by (detC/)/j = 
QiiQjj — Gij- This tensor appears in intermediate, gauge- 
dependent steps in our computations, but not in the final, 
gauge-invariant results. 

The first term proportional to the electric field was first 
derived in the original paper of Karplus and Luttingei'^ on 
the anomalous Hall effect in ferromagnets. 
The appearance of interband couplings at the second-order 
gradient expansion is not inconsistent with the quantum 
adiabatic theoremP^ The latter allows transitions to other 
states to occur with a probability quadratic in the adi- 
abatic parameter, which are space-time gradients in our 
problem. See appendix [I] 

One should distinguish this from the case when Ap — 
but the texture is topologically trivial. 
''^ Furthemore, the Bianchi identity which was used in the 
derivation of Eq. ( 19 1 will acquire source terms at these 



points, cf. appendix [E] 

In the wave-packet formalism of Ref.[71 the Hamiltonian 
is expanded to linear order in gradients, and therefore to 
linear order in electromagnetic fields. 

Strictly speaking, the ground state is a Fermi sphere in 
the p space, and the corresponding distribution is slightly 
shifted by a in the k space. 



Left-handed (LH) Weyl fermions [Eq. ( 49 1 with a minus 



sign] has a zero mode with the opposite chirality so that 
particles leave the system under spectral flow. Therefore, 
Weyl fermions may occur in pairs with opposite chiral- 
ity, so that particles are transferred between the crossing 
points in momentum space, conserving particle number. 
This creates a nonequilibrium distribution in momentum 
space, which would relax back to equilibrium through scat- 
tering between crossing points.'^ 

At this order, second-order derivatives of the diagonal dis- 
tribution function will appear, and the transport equation 
cannot be interpreted in terms of continuous quasiparticle 
motion in phase space. 



The off-diagonal, second order term in Eq. (Ill are neces- 
sary only if we want to solve for n'^' . 



